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Abstract: Graph coloring is used in wireless networks to optimize network 
resources: bandwidth and energy. Nodes access the medium according to their 
color. It is the responsibility of the coloring algorithm to ensure that inter- 
fering nodes do not have the same color. In this research report, we focus on 
wireless sensor networks with grid topologies. How does a coloring algorithm 
take advantage of the regularity of grid topology to provide an optimal periodic 
coloring, that is a coloring with the minimum number of colors? We propose 
the Vector-Based Coloring Method, denoted VCM, a new method that is able 
to provide an optimal periodic coloring for any radio transmission range and for 
any fa-hop coloring, h > 1. This method consists in determining at which grid 
nodes a color can be reproduced without creating interferences between these 
nodes while minimizing the number of colors used. We compare the number of 
colors provided by VCM with the number of colors obtained by a distributed 
coloring algorithm with line and column priority assignments. We also provide 
bounds on the number of colors of optimal general colorings of the infinite grid, 
and show that periodic colorings (and thus VCM) are asymptotically optimal. 
Finally, we discuss the applicability of this method to a real wireless network. 

Key-words: Graph coloring, wireless sensor networks, grid, periodic coloring, 
Vector-based Method, valid coloring, optimal coloring, number of colors bounds. 
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Coloriage des reseaux de capteurs en grille: La 
methode de coloriage basee sur les vecteurs 

Resume : Le coloriage des graphes est utilise dans les reseaux sans fil afin 
d'optimiser les ressources du reseau: la bande passante et l'energie. Les noeuds 
du reseau accedent au medium en fonction de leur couleur. L'algorithme de 
coloriage a la charge d'assurer que les noeuds intereferents n'aient pas la meme 
couleur. Dans ce rapport de recherche, nous nous concentrons sur les reseaux 
de capteurs sans fil ayant une topologie en grille. Comment un algorithme 
de coloriage peut profiter de la regularite de cette topologie pour realiser un 
coloriage periodique optimal, c'est a dire un coloriage avec le nombre minimal 
de couleurs? Nous proposons la Methode des Vecteurs, notee VCM, une nouvelle 
methode qui est capable de fournir un coloriage a h sauts, periodique et optimal 
quelque soit h > 1 et quelque soit la portee radio. Cette methode consiste a 
determiner quels noeuds de la grille peuvent utiliser la meme couleur sans creer 
des interferences entre ces noeuds tout en minimisant le nombre de couleurs 
utilisees. Nous comparons le nombre de couleurs utilisees par la Methode des 
Vecteurs, a celui obtenu par un algorithme distribue qui affecte les couleurs selon 
la priorite donnee par la ligne ou la colonne. Nous fournissons aussi des bornes 
sur le nombre de couleurs des coloriages generaux optimaux de la grille, et 
prouvons que les coloriages periodiques (et done VCM) sont asymptotiquement 
optimaux. Enfin, nous discutons 1' applicability de cette methode a un reseau 
sans fil reel. 

Mots-cles : Coloriage des graphes, reseaux de capteurs sans fil, grilles, colo- 
riage periodique, Methode des Vecteurs, coloriage valide, coloriage optimal, 
bornes sur le nombre de couleurs. 
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1 Context and motivation 

Node coloring consists in assigning colors to nodes using a minimum number of 
colors such that two interfering nodes do not have the same color. It has been 
proved that the 1-hop coloring problem in pQ and more generally the h-hop 
node coloring problem with h an integer > 1 is NP-complete [UGH- That is why, 
approximation algorithms are adopted. Among them, FirstFit [7] assigns to the 
uncolored node with the highest priority the first available color. Unfortunately, 
the number of colors obtained depends on the node coloring order given by the 
priority. A whole class of results expresses properties related to the worst-case 
performance of approximation algorithms: they typically prove that, for any 
input graph of a given family, the coloring obtained by a given algorithm uses 
at most a times the optimal number of colors. Such an algorithm is denoted an 
a-approximation algorithm. 

Many coloring algorithms have been designed to be used in wireless net- 
works to make the medium access more efficient. Among them, we can cite 
TRAMA [5], FLAMA 0, ZMAC [TO] which is based on the DRAND [TI] col- 
oring algorithm, TDM A- ASAP [TO], FlexiTP [TO] and SERENA [TO]. In short, 
time slots are assigned per node color: two nodes with the same color can 
transmit in the same time slot without interfering. 

In this research report, we focus on wireless sensor networks with grid topolo- 
gies. These topologies are regular. Wireless sensors organized according to a 
grid topology are a natural choice and exist in many real applications. Often, 
the grid organization is easy to deploy, and is efficient in terms of coverage, 
connectivity and management. For instance, a grid topology is one of the best 
methods to ensure sensor coverage for surveillance [HE]. It also helps to collect 
measurements with a uniform spatial sampling such as for instance in precision 
agriculture and irrigation as in [BJ: when physical phenomena are numerically 
modeled, measurements from a grid pattern may be a direct input or directly 
compared to the output of equations solved with the finite element method on 
a grid — without requiring additional sensor localization and numerical mea- 
surement interpolation. 

Given such grid topologies, the question (and our problem statement), is: what 
is the method to color grids while maximizing the color reuse? Naturally, a first 
step towards answering the question, is to apply existing algorithms to a sample 
of grids and analyze the output for some sample topologies. 
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Table 1: Number of colors for various transmission ranges, grid sizes and priority 
assignments., for 3-hop coloring 



Grid size 


priority assignment 


colors 
for range = 1 


colors 
for range = 2 


10x10 


line or column 


8* 


30 


diagonal 


8* 


28 


distance to origin 


8* 


30 


random 


13 


36 


20x20 


line or column 


15 


33 


diagonal 


8* 


29 


distance to origin 


8* 


30 


random 


14 


41 



We have performed simulations of a 3-hop coloring, because 3-hop coloring is 
needed when broadcast and immediate acknowledgement of unicast transmission 
are required. The algorithm colors nodes according to their priority order among 
their neighbors up to 3-hop. Node priority is given by one of the following 
heuristics: the position of the node in the line, column, diagonal, its distance to 
the grid center or random. Results are summarized in table Q] The '*' symbol 
highlights the optimality of the number of colors used. The reader can refer 
to [3] for the proof of the optimal number of colors needed to color some grids 
with various transmission ranges. Different radio ranges are simulated. Mostly, 
we observe that no priority assignment tested provides the optimal number of 
colors for any grid size and any radio range. 

This prompts the following questions: Does a priority assignment exist in grid 
topologies such that the number of colors does not depend on node number but 
only on radio range? Is it possible to take advantage of the regularity of grid 
to design a /i-hop coloring algorithm able to find the optimal number of colors 
for any grid with any transmission range value? Moreover, does a periodic color 
pattern that can tile the whole topology for a given radio range, exist? What 
is the optimal number of colors? Is the optimal number of color reached by 
periodic colorings? 

In this research report, we answer these questions by proposing a coloring 
method, the Vector-Based Coloring Method, denoted VCM which provides an 
optimal and valid coloring of grids. 

The research report is organized as follows. In Section [5J we introduce 
some definitions and notations. We present the problem statement and give an 
intuitive idea of VCM and introduce its main components. In Section [3J we 
focus on periodic coloring. In Section El we consider periodic colorings where 
any color is periodically reproduced, and show how VCM computes the color of 
an node in the grid. In Section [SJ we restrict or study of periodic coloring to 
valid ones. We propose two methods to check the validity of periodic coloring. 
Section [5] allows us to characterize an optimal coloring and give upper and 
lower bounds on the number of colors, and shows that periodic colorings (and 
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therefore VCM) are asymptotically optimal when the radio range increases^. We 
then show how to reduce the research space of the generator vectors produced 
by VCM in order to reduce the complexity of VCM that we evaluate (which is 
polynomial time). We conclude this section by giving the method to find the 
optimal valid generator vectors. In Section we summarize by giving rules to 
apply VCM in practice. Results obtained by VCM are given in Section [5] and 
compared with those obtained by different heuristics. Finally, in Section we 
discuss how to apply VCM in real wireless sensor networks. We conclude in 
Section 1101 In the Annex, the reader will find some mathematical results and 
grid properties useful for grid coloring. 

2 Problem Statement and Overview 

2.1 Notations, Definitions and Assumptions 

In the remaining of this research report, we use the following notations: 

• UV denotes a vector of extremities the nodes U and V. 

• \UV\ denotes the norm, or length of the vector UV, and d(U,V) the 
euclidian distance between nodes U and V . 

• det(UV\, UV-i) is the determinant of the two vectors. 

• UV\ ■ UV2 denotes the scalar product of the vectors. 

• A(UV\_, UV2) denotes the lattice having as a base the vectors (UVi, UV2). 
We call UV\ and UV2 the basis (or generator) vectors. V(UV\, UV2) is the 
fundamental parallelotope associated to A. In 2 dimensions, V is a par- 
allelogram. Moreover, the number of nodes in V{UV\,UV2) called lattice 
determinant is given by det{UV\,UV2)- These notations and definitions 
are adopted from [2]. We also denote Pu(UVi, UV2) the fundamental 
parallelotope translated at node U. 

• For x,y in Z 2 , with y 7^ 0, we denote "x modulo y" the integer z in 
{0, 1, . . . \y\ — 1} such that we have x = z(mod\y\). 

Moreover, the following definitions and assumptions are used: 

• Nodes: The nodes are disposed in a grid. Without loss of generality, we 
assume that the grid step is 1, hence the set of nodes is identified by 1? . 

• Neighborhood: R denotes the radio range and is a real number > 1. 
The underlying assumption is a unit disk model, assuming that nodes U 
and V can communicate if their euclidian distance is lower than R and 
are able to communicate directly in both directions. Hence, the set of 
neighbors of a node U is: 

J\f(U) = {V e Z 2 I < d(U, V) < R}. 

1 showing that specifically for grids, there are better techniques than already known general 
a-approximations (even for unit disk graphs, such as 1191 ). since here a = 1 + 0(1/ R) when 
R — >• 00 
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• Coloring: A coloring </> is a mapping from the nodes Z 2 to a set of colors, 
identified with the set of integers {0, . . . k — 1}, with k a positive integer. 

• Valid Coloring: a coloring is said to be a valid h-hop coloring, with h 
an integer > 1, when all nodes that are less or equal to h-hop away are 
assigned different colors. 

• Periodic Coloring: A coloring is denoted "periodic" if the set of nodes 
with the same given color, is the lattice A(vi,V2) for some fixed vectors 
V\, V2 of Z 2 , after it is translated by some amount. In this paper, we 
consider a "strict" definition of periodicity because colors might otherwise 
be periodically repeated with a pattern generating more than a single 
lattice for instance. 



2.2 Problem Statement 

Our general goal is the following: Find a valid h-hop coloring of the nodes of 
the infinite grid Z 2 , with a minimal number of colors for h > 1. 
Because the set of colorings of Z 2 is infinite, in this paper we restrict ourselves to 
colorings exhibiting a periodic color pattern (following the definition of periodic 
coloring of section 12. ip . 

Problem Statement: Find one of the valid periodic h-hop color- 
ings with the minimum number of colors for h > 1. 



2.3 Overview 

In this paper, we propose the periodic Vector-Based Coloring Method, denoted 
VCM that answers the problem statement in polynomial time. This method 
consists in producing a periodic coloring obtained by tiling a color pattern. This 
color pattern is generated by two vectors verifying some conditions that will be 
detailed in the paper. 



As illustrated in Figure [TJ VCM is composed of the following components: 



Parameters R,h VCM 

I Optimal Vector Search 

1 ► (VCM-OVS) 



Parameters R,h 

Generator ». 

vectors u ,u 



Optimal 
■ generator . 
vectors 



VCM 
Validity Check 
(VCM-VC) 



Valid: 
yes or no 



Method 1 



Method 2 



Generator 
vectors u ,u 

1 2 

i 

VCM 
Node Color 
Computation 
(VCM-NCC) 

~ *~ 
Colors of nodes 

o o • 



Figure 1: Components of VCM 
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CI. VCM-NCC A periodic coloring based on two generator vectors. This 
component determines how any sensor node selects its color. 
C2. VCM-VC An algorithm to check the validity of the periodic coloring: the 
VCM incorporates two methods to verify the validity of the coloring for a given 
couple of generator vectors. 

C3. VCM-OVS An algorithm to search for optimal generator vectors, that is, 
yielding the minimum number of colors. Indeed, we limit the set of candidate 
vectors to find the vectors providing the optimal number of colors. 



2.4 Intuitive Idea 

The intuitive idea of the method is as follows. As the grid topology presents a 
regularity in terms of node position, VCM produces a similar regularity in terms 
of colors and generates a color pattern that can be periodically reproduced to 
color the whole grid. Our aim is to find such a color pattern that minimizes the 
number of colors used. More precisely, given a colored node U, we determine 
where its color can be reproduced. The method starts by the choice of two 
vectors UV\ and UV2 such that V\ and V2 use the same color as U. Of course, 
the vectors UV\, UV2 must provide a valid h-hop coloring. The color pattern is 
given by the set of colors in the finite parallelogram V(UVi, UV2) translated at 
U. Hence, the color of U is repeated also at any node W where UW is a linear 
combination of UV\ and UV2- 

In the Sections from @] to [51 we detail the components of VCM. 



3 Periodic Coloring 

Figure [3 presents an example of a periodic 3-hop coloring of a 10 x 10 grid with 
a transmission range R = 1. 

o • • tp,t • • O • 

o o • ®' o \m 000 

• • • • » • • 

•©©^•0 0^#O 

o • • •» • • p • 

o o m@'o c\# 0' 

• •v i gt'# • • V v '« • • 

• O O TS • Q^OO • O 
O • • • • • • 

OO#O00#O00 

Figure 2: Example of periodic 3-hop coloring (R=l) 
The principles of the periodic coloring are: 
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PI. (Generator vectors) The two vectors UV\ and UV2, if linearly indepen- 
dent, generate the parallelogram V(UVi, UV2) of the color pattern. 
P2. (Parallelotope color unicity) Inside T , (UVi,UV2), there is no color 
reuse. 

P3. (Lattice color repetition) Because the periodic coloring is obtained by 
repeating the color pattern, the nodes W such that UW is a linear combination 
of UVi and UV2 have the same color as U. Precisely, the set of nodes W having 
the same color as U forms a lattice of Z 2 with generator vectors UVi and UV2: 
the vector UW can be written as UW = aUV\ + PUV2 with a and /3 in 1? . 
P4. (Coordinate-based color computation) The grid can be seen as a 
tiling with the pattern VuiUVi, UV2). Thus, each node W belongs to its own 
parallelogram, and has coordinates relative to this parallelogram. Consequently, 
W has the same color as any node having the same coordinates in its own par- 
allelogram. 

To be applied to a wireless sensor network, these principles have to be en- 
hanced. For instance, since nodes having the same color can simultaneously 
access the wireless medium, validity of the coloring must be verified. Moreover, 
to ensure an efficient usage of the bandwidth, the number of colors used should 
be minimized. These criteria are taken into account in our work and progres- 
sively introduced in the paper. 

In the following we detail the components of VCM. As previously said, we 
only consider grid colorings that periodically reproduce a color pattern. 

4 VCM: Node Color Computation (NCC) 
4.1 Assigning Colors to Nodes 

The node color computation (NCC) component of the VCM method takes as 
parameters two generator vectors U\,u-x (computed as in Section[5|). Let (x\, y±) 
and (#2,2/2) be their coordinates and let d = det(ui,U2). Here we define two 
methods to compute the colors. 

Actual computation on an example is provided in Section 14.31 
4.1.1 Method NCC1 

Method VCM-NCC 1: VCM assigns the color of a point W based on its 
coordinates w = (x,y) by computing first the integer quantities ci(w), 02(10) as 
in System (TTJ), 

f c\(w) = (xy 2 - yx 2 ) modulo d , . 

\ Ciiw) = (—xyi +y x i) modulo d 

and then using a bijective mapping between the couple (01,02) and a color 

e{o,i,2,...|d|-i}. 



Remark 1 Defining C2 as c^{w) = (xyi—yxi) modulo d, instead of the defini- 
tion in (Tip, would be more symmetrical compared to C\, and would yield identical 
results (as it is a trivial transformation: 02(11)) — > (—02(10)) modulo d). 
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Remark 2 In the remaining of this report, we will assume that d > without 
loss of generality: indeed, if d < 0, it is sufficient to use the vectors (— Ui,v>2) 
instead of (1x1,1x2) and the results are similar ; notice that c\(w, 1x1,1x2) = 
C\(—W,— 1x1,1x2), etc. (hence change of sign of u\ is equivalent to an origin 
symmetry of the coloring). This avoids minor technicalities on the definition of 
the modulo, integer part, fractional part, when numbers are negative (for which 
definitions are not universal) . 

Property 1 With the previous coloring VCM-NCC1, it is indeed possible to 
define a bijection from (ci(w), c 2 (w)) to {0, 1,2,... \d\ — 1}. Moreover, the col- 
oring verifies principles P3 and P4 defined in Section^ 

Proof: We assume that d > without loss of generality (see Remark [5]) . We need 
to prove that the set {(ci(w),C2(w)) \ w £ Z 2 } has cardinality d and that the 
set of nodes with the same color is exactly the lattice A(ui, u 2 ) translated at w. 

Let W be a grid point of coordinates w = (x,y). Performing a change of 
vector basis in IR 2 from {(1,0), (0, 1)} to {1x1,1x2}, the new coordinates (a, (3) £ 
R 2 of W in A(ui, 1x2) verify w = au\ + fiui and 

dct(ui, U2) 
a = 

_ dct(ui, w) ^ 

v d 

with d = det(ix\, u%). 

Let a' and ft' be the integer parts of a, /3, i.e., a' = [a\ and j3' = \_/3\. 
For arbitrary nonzero integers A, fi, with also fi > 0, we have the identity: 
— = I —J + Amodul °^ . Thus 121 becomes: 



, detiw, U2) modulo d . cUw) 
a = a H = a H — 

, det(ui,u;) modulo d _ , c 2 (w) 

P P ~j ' P 1 ~j 



(3) 



Let W the point with coordinates w' = (a',/3'). W is on the lattice since 
a', j3' are integers, and observe that W is in fact inside the parallelogram of the 
lattice A(wi, M2) placed at node W' (i.e. inside the parallelogram defined by the 
4 points of the lattice: w' , w' + w' +112, w' + U1+U2). Then ([3]) means simply 
that ( Cl y^ , C2 y"^ ) are the coordinates of W relative to this parallelogram (with 
the basis vectors Ui,ti2) . 

Since there is a bijection between the set of coordinates of nodes in a par- 
allelogram of A(wi,ii2) an d the nodes themselves; and since ( cl y"^ , T ) are 
these coordinates, we have the two properties: 1) there are exactly d possible 
values of (ci, C2) (because there are exactly d nodes in the parallelogram), and 
2) no two nodes inside the parallelogram have the same values cl, c2 since these 
are their coordinates, relative to one vertex of the parallelogram. ■ 

Lemma 1 With the color computation given by the System [TJ the color of the 
node U is repeated at the nodes W with coordinates verifying: w = au\ + (3u2, 
for some (a, f3) € Z 2 . 

Proof: Actually, by construction, the color of a node is given by its coordinates 
relative to the parallelogram it belongs to. Hence, the color of any node U is 
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reused at nodes W, such that UW = a,u\ + j3u2 for all (a, /3) € Z 2 , which have 
the same relative coordinates. ■ 
We deduce that VCM-CC1 provides a coloring that is really consistent with 
the principle of the method as described in Section [3] 

4.1.2 Example of bijection for NCC1 

One of the steps of the Method NCC 1, is that a bijection needs to be estab- 
lished between the set of values {(a(w), C2(w)) | meZ 2 } and the set of colors 
{0, 1, . . . , d — 1}: an example of bijection is provided in this section. 

A bijection can be constructed by computing the values of (c±(w), C2{w)) for 
any node in V(u\,U2) in a list, sorting the list by lexicographical order, and, 
finally, setting the color associated with a couple (ci(w), C2(w)) to be its index 
in the sorted list minus 1. 

Example: if (0, 0) appears as the I s * item of the sorted list [as it can be 
proved it will], the color assigned to that couple is 0. Then, for instance, for 
the point W of coordinates w = (xi,yi), in other terms w = u\, we have 
(ci(w), c 2 (ui)) = (0,0) and therefore the color assigned to W is 0. 

Note that, from a pure implementation point of view, it may be difficult 
to enumerate exactly the points of T^iti, u 2 ), but then, instead, it is sufficient 
to enumerate all the nodes in a superset, the bounding box of V(u\,U2), itself 
computed from its four vertices O = (0,0), O + Ui, O + U2 and O + u\ + U2- 
Computing the set of values (ci(w), C2{w)) for the points in the bounding box, 
will yield all possible values for any aieZ 2 . 

4.1.3 Method NCC2 

Method 2 is derived from the first method ; the difference is that it also estab- 
lishes a direct bijection (and hence avoids the need for constructing a bijection 
as in section I4.1.2p : it proceeds to a direct computation of the color based on 
the node coordinates without computing the colors of the other nodes on the 
parallelogram. 

As for Method 1, we assume that we are given two generator vectors Wi,it2 
with coordinates {x\,y\) and (2:2,2/2), anc ^ that we are computing the color of 
node of coordinates w — (x, y). 

In addition to the notations and definitions used previously, we introduce 
the following ones: 

• gi is the greatest common divisor of (xi,yi); g\ = gcd{x\,yx). Similarly, 
92 = gcd(x2, 2/2), with the convention gcd{a 1 0) = gcd(Q, a) = a 

• vi is the vector = -^iti and V2 is the vector = ^7 U 2- Notice that the 
coordinates of v\ and V2 are coprimes. 

• Let d' = det(vi 1 v 2 )- We have: d! = 

• Let (a'(w), P'{w)) £ M. 2 be the coordinates of the node w relative to the 
basis (ui,i>2), when performing a change of basis in M 2 from (1,0), (0, 1) 
to (vi,v 2 ). We have: a'(w) = and p'( w ) = det(^,w) 

• Let x'(w) = [a'(w)J modulo g\ 
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• Let y'{w) = \J}'(w)\ modulo 52 

• Let d(w) = det(w,V2) modulo d' 



Method VCM-NCC 2: Using the notations mentioned above, a color of a 
node with coordinates w = (x, y) is equal to 

c( w ) = c'(w) + d'x'(w) + d'giy'(w) (4) 

which is an integer in {0, 1,2,... \d\ — 1} 

As in remark [2j we now assume that d > without loss of generality. 
The idea of the method VCM-NCC 2 is as follows: 

• In method VCM-NCC 1, we did not identify nor did make use of the 
special structure of E c = {(ci(u>), C2(w)) \ w <E Z 2 }, which is in fact a 
subgroup of (J|) 2 ; hence we did not provide an explicit mapping from 
E c to {0, 1, 2, ... d — 1}. This is the problem that method 2 addresses: 

• In method VCM-NCC 2, by dividing the vectors ui and 1x2 by the re- 
spective gcd of their coordinates, we obtain vectors vi,t>2, and we find 
that applying method VCM-NCC 1 with these vectors, the set of val- 
ues {ci(u>, v\, V2), 02(11}, vi, V2) I w £ Z 2 } has good properties, allowing the 
construction of a direct bijection (see lemma [5]). 

• However using method VCM-NCC 1 with vectors vi,V2 yields a coloring 
of Z 2 repeated by translation by V\ and «2 ; whereas we wanted a coloring 
repeated by translation by u\ and U2 (which are larger). Notice that the 
problem here is that method VCC-NCC 1 with generator vectors v\ , V2 is 
quite possibly an invalid /i-hop coloring of Z 2 , even if u\,U2 are vectors 
defining a valid h-hop coloring. For this reason, method VCM-NCC 2 
constructs a coloring based on method VCM-NCC 1, but modified, by 
tiling the parallelogram T , (v\,V2) several times inside the parallelogram 
V (1*1,1*2) and changing the colors in each internal tile (see figure [3]). 

Method 2 is in fact the combination of two functions; c(w) = 63(03(10)) with: 

• The first one, c 3 , fromZ 2 to {0,1, ..<f - 1}x{0,1, ..g x - l}x{0, 1, ..g 2 - 1}, 
which associates to w the value 03(10) = (d (w) , x' (w) , y' (w)) . 

• The second is bijection 03 from {0, 1, . . . d' — l}x{0, 1, . . .g\ — l}x{0, 1, . . .32 
to {1, 2, ... d} by transforming (j, k, I) — > 63 (j, k, I) = j + d'k + d! g\l ; note 
that d = g\g2d' by definition. 

The transformation 63 is obviously a bijection from its domain to its codomair|^|, 
and therefore the crux of method VCM-NCC 2 is the function C3. 

We start with the following lemma for the function d(w), which is in fact 
ci(w, vi,v 2 ): 

2 informally, it is similar to a transforming the time of the day expressed as h:m:s - h 
hour m minutes s seconds - into the time of the day expressed as seconds from midnight as 
(h, m, s) -> 3600 xh + 60xm + s 
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Figure 3: Method VCM-NCC 2 



Lemma 2 Let C\{w, V\, 1)2) be the value "c\(w) " obtained when applying method 
VCM-NCC 1 with generator vectors (1)1,1)2), and define C2{w,V\,v%) similarly. 
Then the function w — > c\(w, V\, V2) is a direct coloring of the nodes w, where 
the color is in {0, 1, . . . d! — 1}. 

This lemma establishes that for the specific generator vectors (v 1,1)2), c\ com- 
putes directly a color and C2 is essentially redundant. 
Proof: We prove this lemma in two steps: (SI) proves that the set {ci(u>, v\, V2) \ w€l?} 
covers exactly integers in {0, 1, . . . d' — 1} , and (S2) the coloring provided by c\ 
has the same properties as VCM—NCCl. To prove this, we prove that there is a 
bijection between {ci (w, v\, V2) \ w £ Z 2 } and E c = {(ci(w, 1)1,1)2), C2(w, v\, 1)2)) \ w € Z 2 }. 

(SI) Consider the vector 112- By definition, 172 is the gcd of its coordinates 
(^2,2/2)- Applying Bezout's identity ([20]), there exist integer^ (A,/i) £ Z 2 , 
such that: 

Xx 2 + /iJ/2 = 92 (5) 

Let wo be the vector of coordinates (/i, — A) and let (x' 2l y'2) be the coordinates 
of V2 (recall that i>2 is defined as 1)2 = ~^ u 2)- Because V2 = j^ u 2, we have 
and y! 2 



■f 2 

92 



U and hence equation (JSJ) becomes: 



Ax 2 + iiy'2 = 1 



det(w ,v 2 ) 
=> c 1 (w ,v 1 ,v 2 ) 
Since c\ is a linear map: cx(2iuo, v\, V2) - 



1 



2, etc. and by applying it to 



0, wq, 2u>o, ...(df — l)wo, we get exactly: 

{ C1 (W,V U V2) \ wEl 2 }= {0,1, ...d! - 1} 



'this is even true when £2 = or yi = with our convention gcd(a, 0) = gcd(0, a) = a 
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(S2) Let E c be the set E c = {(ci(w, v\, V2), Oz{w, v\, v%)) \ w € Z 2 }. From 
properties previously proven for method VCM-NCC 1 (including property [1} , 
we know that there exists a bijection from E c to the set {0, 1, 2, . . . d' — 1} ; in 
other words, that its cardinality verifies \E C \ = d' . 

Now consider the projection p x : (x, y) — > x. The previous relation shows 
that p x {E c ) = {0, 1, . . .d' — 1} hence it is a bijection because it is surjective 
between two sets of cardinality d' . 

Therefore we have explicitly found one bijection (precisely: k,l — > p x (k,l)) 
compared to the method VCM-NCC 1 where we only established the existence 
of such a bijection. The end result, is that when applying VCM-NCC 1 with 
the generator vectors (v\, V2) and with this bijection p x , the color of a node w 
is Ci(w, vi, V2). The coloring provided by c\ inherits from VCM-NCC 1, all the 
properties listed in Section [3] I 

The next step is to build a coloring repeated by translation of u\ and U2, from the 
coloring VCM-NCC 1. As illustrated in Figure [31 the idea is that V\,V2 are the 
basis of the parallelogram V{v\, V2) that can be tiled to cover the parallelogram 
T D {ui,U2) ; this is possible because u\ = g\V\ and 112 = 32^2 (where g\ and 
<?2 are integers). Inside a parallelogram V(vi,V2), the nodes are assigned the 
same color as given by VCM-NCC 1 with vi, V2, except that an offset is added, 
depending in which tile the node is located. 

Concretely, in Figure [3J the parallelogram V Ul (ui,U2) at the right side of 
the picture is tiled with 6 smaller parallelograms which are versions of V(vi, V2), 
filled by different shades of gray. The nodes in one small gray parallelogram are 
allocated the colors 0, 1, 2 . . . d'— 1 ; the ones of the next small gray parallelogram 
are allocated the colors d' ,d' + 1, . . . 2d' — 1; and generally the nodes of the k th 
parallelogram are allocated the colors kd' , kd' + 1, . . . (k + l)d' — 1. The actual 
coloring is then something like the colored nodes at the left of the Figure [3J 

Thus, the central idea is to identify in which sub-parallelogram a node is 
located: this is given by x' (w) , y' (w) , which are in fact the coordinates of the 
smaller parallelogram on the basis (t>i,i>2) counted relative to the larger paral- 
lelogram. For instance, in Figure [3J the point w is inside a sub-parallelogram of 
coordinates (x'(w), y'(w)) — (1,2) relative to the bigger parallelogram. 

This is the informal idea behind the quantities defined previously: in the fol- 
lowing lemma, we now prove formally that VCM-NCC 2 is a coloring satisfying 
the desired properties. 

Lemma 3 The coloring provided by the method VCM-NCC 2 satisfies the prin- 
ciples PI, P2, P3 and P4 from Section^ 

Proof: 

Almost all the properties can be derived from identifying the set of nodes 
with identical colors. 

Let w\ and W2 be two nodes with the same color, hence, with c(w\) = c(u>2). 
Because c = 63 o C3 and 63 is a bijection, this implies that 

c'(wi) = c'(w 2 ) and x' {wi) = x'(w 2 ) and y'(wi) = y'(w 2 ) (6) 

The part c'{wi) = c'{vj2) implies that w\ and W2 have the same color in the 
method VCM-NCC 1 applied with vectors v\ , «2 , therefore, by properties of this 
method: 

W\ — W2 = A«i + /ii>2 for some (A, n) € 1? 
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Let us decompose A as quotient and remainder modulo g\ (resp. fi, mod- 
ulo g 2 ): A = #1 gi + n and \i = g 2 q 2 + r 2 where r x G {0, 1, . . .51 - 1}, r 2 € 
{0,l,...fl2-l} and ( qi ,q 2 )eZ 2 . 

The previous equality becomes: 

W\ - w 2 = {giqi + ri)vi + (g 2 q 2 + r 2 )v 2 (7) 
Now developing the part x'{w\) = x'(w 2 ) from ([5]), we get: 
|_c/(wi)J modulo g\ = \ ol (w 2 )\ modulo g\ 

=>■ [a'(wi) — c<'{w 2 )\ modulo gi = or [a'(w 2 ) — a'(wi)\ modulo g\ = 

Which equality is true depends on which one of a'(wi) and a'(w 2 ) has the 
largest fractional part. Assume it is a'(wi), then we have: 

[a'(wi) — a'{w 2 )\ modulo g± = 
y det( Wl -w 2 ,v 2 ) j modmQ ffi = o ( by dcnnition of a ') 

=^ L(.9i<?i + ri) deti - v ^' v ^ j modulo gi = (using cqQ 
==> ri = (using the definition of d') 

In the same way, developing y'(w\) = y'(w 2 ) from ([5]), yields r 2 = 0. Since 
Mi = giV\, u 2 = g 2 v 2 , we can now rewrite Eq ([7]) as: 

wi — W2 = qiU\ + (72^2 for some (qi, g 2 ) G Z 2 

The opposite is true: if w± — w 2 verify this equality. ui\ = w 2 + qiUi + q 2 u 2 , 
and then one can easily check that c'(w\) = c'(w 2 ) and similarly: 

x(wi) = L a '( u '2 + qi u i + 92^2)] modulo gi 

det(w 2 + qiui + q 2 u 2 ,v 2 ) 

= [ J77 ^ J modulo gi 

aet(vi, v 2 ) 

= \_a (w 2 ) + 9i<?iJ modulo gi 

= [a'(w 2 )\ modulo g\ = x'(w 2 ) 

Likewise y'{wi) = y'(w 2 ) and thus finally c(wi) — c(w 2 ), that is, the two nodes 
of coordinates w\ and w 2 have the same color. 

As a result, we have established that the nodes with identical colors are 
exactly the points located on a lattice generated by vectors v\ , v 2 . This is 
exactly the principle P3, and actually implies the principles PI, P2 and P4. 



Remark 3 In this section, we selected c'(w) as c'(w) = C\{w,Vi,v 2 ). Alterna- 
tively, one could select c'(w) = c 2 (w,vi,v 2 ). Notice also that once the choice is 
made, lemma\^does not require that coordinates are coprimes for both vectors V\ 
and v 2 . Indeed, for the choice d = c\, ( respectively c' — c 2 ) only v 2 ( respectively 
v\) is required to have coordinates that are coprime. 
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4.2 Computing the Number of Colors 

The number of colors used in a periodic /i-hop coloring is given by the next 
Property. 

Property 2 For any node U , the color pattern defined by the two generator 
vectors ui andu 2 meeting the aforementioned principles contains exactly \x1y2 — 
X2yi\ colors where (xi,yi) and (x 2 ,y 2 ) ore the coordinates of u\ andv.2- 

Proof: By definition, no two nodes within the parallelogram defined by u\ and u 2 
use the same color. Hence the number of colors is equal to the number of nodes 
in this parallelogram. Moreover, as we said, the number of nodes in V(ux,U2), 
called lattice determinant, is equal to the absolute value of det(iti,it2) |14| . 
Hence the property. ■ 



4.3 Example of Color Calculation 

In this section, we illustrate the color calculation, using the example of Figure [3] 
In Figure [3l we have the following coordinates for Ui, 112, w: 

• ui = (6,2) and u 2 = (—3,6) 

• w = (8,9) 

Applying VCM-NCC 1 with the generator vectors ui,u 2 , we get: 

• Number of colors = 42 

• (ci(w, ui,w 2 ),c 2 (w, ui,u 2 )) = (33,38) 

• Using the example bijection of Section l4.1.21 ci(w), c 2 (u;) is the 36 th value 
in the sorted list of possible values, hence color(w) = 35 

Applying VCM-NCC 2 with the generator vectors ui,u 2 , we get: 

• Number of colors = 42 

• 5i = gce?(coords of ui) = 2 and then V\ = (3, 1) 

• .92 = gcd(coords of u 2 ) = 3 and then V2 = (—1, 2) 

• d{w) = 4, x'{w) = 1, y'(w) = 2 

• Color of w. c(w) = 39 

For reference, the colors are computed internally from the coordinates w on 
different basis: w = (1 + §§)«! + §§«2 and w = (3 + |)ui + (2 + | )v 2 



5 VCM: Validity Check (VC) 

As defined previously, a /i-hop coloring algorithm is valid if and only if no two 
nodes that are at less or equal to /i-hop from each other use the same color. 

The node color computation algorithm of VCM (described in Section^]) takes 
as input two generator vectors U\, and U2, and gives the color of each node. In 
this section, we will assume that such two vectors are given and fixed, and we 
present two methods for checking beforehand whether the coloring induced by 
these vectors is a valid coloring. 
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5.1 Method VC1: Verification around Origin 

Method VC1: For each node W in the /i-hop neighborhood of the origin node 
O, we compute the color of this node based on the given generator vectors i*i 
and U2- If W has the same color as O, then we conclude that the vectors u\ and 
1*2 do not provide a valid coloring. Otherwise, if the color of O, is not repeated 
at any point W in its /i-hop neighborhood, then the coloring is valid. 

The idea of Method @] is based on the following fact, proven in this section: if 
there is a color conflict between any two nodes V\ and Vi in A(ui,w 2 ) (Vi and 
V2 have the same color despite they are at less than or equal to h hops), there 
will be a color conflict in the /i-hop neighborhood of the origin O. 
We set d = det(ui, 112)- 

Lemma 4 If two nodes Vi and V2 with coordinates V\ , V2 in Z 2 have the same 
color, then the color of the origin node is repeated at the node W of coordinates 
v 2 -v 1 . 

Proof: The functions ci , C2 computed from System [1] are actually linear modulo 
d. That is, if W is the node with coordinates v\ — V2, and w is the vector of 
nodes extremities the origin and W, we get: c\(W) = c\ (Vi) — 01(^2) modulo d. 
Hence, if ci(Vi) = Ci(V2) we have Ci(W) = c\(0). This is true also for C2, hence 
the lemma. ■ 

Theorem 1 If none of the nodes inside the h-hop neighborhood of the origin 
node O = (0, 0) has the same color as O itself, then the coloring is valid. 

Proof: By contradiction: assume that the coloring is invalid, which implies that 
two nodes V\,V2 at less or equal to h hops have the same color. Then from 
Lemma |U the node W such as OW = V1V2 has the same color as O. Notice 
that the distance in terms of hop number between O and W is the same as the 
distance between V\ and V2- Hence we have found a color conflict between O 
and a node W which is at less than h hops from O. Hence the theorem. ■ 
Theorem[T]proves that Method VC1 is a correct method for checking whether 
two generator vectors yield a valid /i-hop coloring. 



5.2 Method VC2: Verification in a Few Points 

Method VCQ] requires 9(i? 2 ) verifications when R — > 00. In the following, we 
propose Method VCI21 usable when R > \pl and requiring only a bounded num- 
ber of verifications. Method VQUperforms a check on a few nodes on the lattice 
K(u\,U2) to guarantee that the /i-hop coloring associated to U\,u<2 is valid. 
This method is based on Gauss lattice reduction [T5] (see the Annex for more 
details): u\ and U2 should be first reduced, and hence verify the Equations (fTUj) . 
Method VC2: The nodes with the same color as the origin are on the lattice 
A(iti,U2): this method verifies that these nodes are at least {h + l)-away from 
O, in which case the coloring is valid. However, not all grid nodes need to be 
checked. It is sufficient to check only nodes W in A(u\,U2) with coordinates 
a, P on the basis {1*1,1/2}, such that |a| and < fx(R), with fi(R) = ^|r^7jy- 
The coloring is valid if and only if these nodes are strictly more than h hops 
from the origin node. 
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This method is based on the following theorem. 

Theorem 2 For R > v2, the coloring provided by two reduced vectors ui,U2 
is valid if and only if: 

for all a, f3 integers verifying \a\ < n(R), and \f3\ < /J,(R), the node with coordi- 
nates (a, f3) on the basis \u\, 1*2} is at strictly more than h hops from the origin 
node O, where = 3 ff^ ■ 

Proof: The property comes from the fact that the points on the lattice A(ui,U2) 
are "far" from the origin node, because the vectors u\,U2 are reduced. 
Indeed, Lemma [§] (see the Annex) means that any node on the lattice with 
coordinates au\ + (3u2, with \a\ or |/3| > n(R) can reuse the color of the origin 
node O because they are at strictly more that h-hop from O (provided that the 
points of coordinates u\ or U2 are themselves strictly more than /i-hop away 
from O). Hence, to check the validity of the coloring provided by VCM, it 
is necessary and sufficient to check that for all \a\ and \f3\ < fJ,(R), nodes of 
coordinates a, f3 in the lattice A(ui,U2) are strictly more than h hops away 
from the origin of the lattice. This check includes checking the validity of u\ 
and U2 themselves (cases (a, (3) = (1,0) and (a, /3) = (0, 1)) 

Notice that for dense grids (i? —> 00), \i — > 1.15 . . .. This small bound reduces 
the set of nodes to be checked in order to verify the validity of the coloring for- 
given vector candidates. 

In fact for R > g 3v ^ , that is for R > 3.3461, we have [i < 2, hence only 4 
points need to be checked (considering symmetries): u\ (with a = 1,(3 = 0), U2 
(with a = 0,/3 = 1), u\ + 112 (with a = 1, (3 = 1), ui — 112 (with a = 1, (3 = — 1). 

However, Method [T] is applicable for any radio range R, whereas Method^] 
requires (R > y/2). 



6 VCM: Optimal Vector Search (OVS) 

To achieve an optimal spatial reuse, the coloring algorithm should minimize the 
number of colors used. For VCM, our aim is to judiciously choose the generator 
vectors u\ and 112 in order to reduce the number of colors used to color a grid. 
However, by default, the infinite lattice 1? is a possible set for candidate vectors. 
So, to find the optimal vectors in a small set, our approach is as the following: 

• We determine the upper and lower bounds on the number of colors needed 
in a h-hop coloring of the grid, for h > 1. 

• Because for any couple of initial generator vectors, reduced vectors always 
exist, it is sufficient to search for some optimal vectors in the space of 
reduced vectors defined by System (TT5]). 

We will show in Section 16.21 how to bound the set of candidate vectors us- 
ing the properties of lattice reduction and the upper and lower bounds on the 
number of colors to decrease the complexity of the search for the optimal vectors. 
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6.1 Bounds on the Number of Colors in Colorings 

In this section, we prove that the number of colors of optimal colorings when 
R — >• oo is shown to be asymptotically ^h 2 R 2 + O(R), from the combination 
of two bounds. 



6.1.1 Lower Bound 

For the lower bound, we have the following theorem that is valid for any coloring, 
not just periodic colorings. It uses known results on circle packings. 

Theorem 3 The number of colors required to color an infinite grid with R > \/2 
is at least ^h 2 (R- V%) 2 ■ 

Proof: Consider h-hop coloring of the grid Z 2 . Consider a fixed color c, and 
now let S c be the set of nodes having this color. 

We first establish a lower bound of the distance of nodes in S c . Let us define 
p = (R — s/2)h. Consider two nodes A, B of S c . By contradiction: if their 
distance verifies d(A, B) < p, from Lemma [51 they would be at most h-hop 
away, contradicting the definition of a ft,-hop coloring. Therefore, all nodes of 
S c are at a distance at least p from each other. 

Now consider the set of circles C of radius ^p and whose centers are the nodes 
of S c . The fact that any two nodes of S c are distant of more p, implies that 
none of the circles in C overlap. Hence C is a circle packing by definition. From 
the Thue-Toth theorem [T71 \TE\ establishing that the hexagonal circle packing 
is the densest packing, with a density of 7^7=, we deduce that C must have a 
lower or equal packing density. This implies an upper bound of the density of 
set S c of centers of the disks of , 1 _- . 

(p/2) 2 v / 12 

Because each color yields a set of nodes with at most this density, it follows 
a lower bound of the number of colors that is the inverse of this quantity, hence 
the theorem. ■ 



6.1.2 Upper Bound 

For an upper bound, when R > v2 , we construct explicitly a periodic coloring; 
more precisely we construct two vectors v\ , V2 yielding a valid coloring with 
VCM. As a result, optimal colorings and optimal periodic colorings must have 
a number of colors which is lower or equal. 

Because the vectors constructed yield lattices which are close to hexagonal 
lattices (when R — > 00), and because it happens that hexagonal lattices yields 
the densest packing (as used in the proof of theorem [3]), the upper bound will 
be somewhat "close" to the previous lower bound. 

We have the following theorem: 

Theorem 4 The number of colors required to color an infinite grid is at most 
^h 2 R 2 + 2hR + (hR + 2)V2. 

Proof: We proceed with a constructive proof, exhibiting two valid vectors 
which yield the result, using an approximation of an hexagonal lattice. 
Figure [4] illustrates how some points V\ and V2 are constructed. 
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Figure 4: Selecting vectors for a near-hexagonal lattice 



• Starting from the point U, the line with an angle ^ with the horizontal 
line is considered, and its intersection with the circle of radius hR yields 
the point B. 

• Next, the closest point of B on the grid having higher coordinates x and 
y than B is sought and is V 2 with coordinates (x 2 ,y 2 ). 

• Then V\ with coordinates [xx,y\) is selected with (xi,yi) = (2x2,0). 

Notice that by construction x\ > hR, and we have a valid choice of vectors 
UVi and UV 2 . 

Denote (7, 8) the coordinates of BV 2 . By construction: < 7 < 1 and 
< S < 1, hence: \BV 2 \ < y/2. Moreover, 

\A%\ = 
< 

Consequently, we can write n c , the number of colors in the associated color- 
ing as: 

n c = det(UV!,UV2) 

= det(UA, UB) + det(AV u UB) + det(U% BV 2 ) 

< det(UA, UB) + \AW\\UB\ + \UV!\\BV 2 \ 

< ^-h 2 R 2 + 2hR+{hR + 2)V2. 



2x 2 - hR, 

2(hRcos(-) + 7) - hR, with 7 < 1 
3 

2 
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Remark 4 An alternate, simpler, choice of vectors is to compute the integer 
A = \hR\ +1, and select the vectors with coordinates u[ = (A, 0) and u' 2 = (0, A). 
The number of colors is higher than for the near-hexagonal previous vectors ( it 
is h 2 R 2 (l + 0(j^))), hence the vectors cannot yield the result of the next section, 
but the vectors can be used for an upper bound when R < \/2. 



6.1.3 Asymptotic Number of Colors 

Theorem 5 The number of colors n c (R) of an optimal periodic h-hop coloring 
for a fixed h verifies: 

n c (R) = ^h 2 R 2 (l + 0(±)) 

when R — > oo. 

Proof: Combining the lower bound and the upper bound of the two theorems [3] 
and @] yields the result. ■ 

Corollary 1 Theorem\3\is true even considering periodic and non periodic col- 
orings. 

Proof: A periodic coloring is a special case of general coloring (including periodic 
and non periodic). Hence, the optimal number of colors in general coloring is 
less than or equal to n c . ■ 

Corollary 2 VCM is asymptotically optimal even considering all possible valid 
colorings (even non-periodic) . In other terms VCM is an (l+g(i?))- approximation 
of the optimal coloring(s) of the grid with some g verifying g(R) — > (precisely: 
g(R) = O(jj)) when R — !> oo. 

Proof: VCM will find vectors with better or equal performance than those in 
Theorem 2J Indeed, in the worst case, these generator vectors for the "near- 
hexagonal" lattice will be selected by VCM. ■ 



6.2 Finding Optimal Vectors 

The complexity of VCM lies in the generator vectors computation and in their 
validity check. In this section, we show how to limit the set of candidate vectors. 
Let u\,U2 be two candidate vectors and be the angle between them. We search 
hmin, hmax (respectively hmim hmax) the lower and upper bounds of the length 
of u\ (respectively U2). 

1. Considering a h—hop coloring, the vectors u\ and u 2 must be valid. According 
to Lemma [51 we have: 

\ui\ > (R- \/2) h and \u 2 \ > {R- V2) h. (8) 

2. As we said, in order to reduce the set of candidate generator vectors, we 
reduce the size of these vectors by using the lattice reduction algorithm of Gauss. 
A consequence of Gauss property [TO] is that: \cos6\ < \, and hence: 

|sin#|>^. (9) 
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3. As shown in Theorem [31 

det{m, u 2 ) <S h = ^h 2 R 2 + 2hR + (2 + hR)V2. 
It results: |iti||it2|| sm9\ < Sh- UsinglHl we get: 
\/3 

— KIM < |ui||w 2 || sin6»| < S h 



(10) 



And as |ui| < Kj, from [TU] we have: 

^\m\ 2 <S h . 

We now separate the cases where R > \/2 and R < \J~2 

6.2.1 Case R > V2 
Using [SJ from [TT] we get: 



111) 



^-\ Ul \(R-V2)h<^\ Ul \ 2 <S h . 



And using [81 in [TOl 



V3, 



'"2 



\(R-V2)h<S h . 



To summarize, the two generator vectors should verify for R > \/2: 

hmm <\Ul\< hmax and hmin < Kl < hmax 



(12) 



(13) 



(14) 



with: 



' l lmin = h(R-V2) 

"lmax A/ J h 

hmin = h(R — V2) 
2 S h 



(15) 



™ V3h(R-V2) 

In practice, to compute the two generator vectors, we determine the upper 
and the lower bounds of the coordinates of u\ , and u-i using the System 1151 
Notice that we can search the vectors in the half plane (y > 0), because if 
U\ and ui are generator vectors, then their symmetric vectors with respect to 
(y = 0) axis are also generator vectors. Consequently, we have: 



hmax — — h 

< yi < h 



max — — "i-rnax 
max 



and 



^2max — — ^2max 
0<V2< hmax 



(16) 
(17) 
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6.2.2 Case R < y/2 

We use li m in = 0. In addition, we can use the bound implied by the vectors 
proposed in Remark @J that is: S s — (hR + l) 2 , which replaces Sh for the 
computation of h max - Then, instead of a fixed bounds for h max and hmax- we 
propose a bound depending on u\, by using (|10j) we have: l2 mal (ui) = 2 ^^ s , 
and because \u2\ > \ui\ we have: h m i n (ui) = \ui\. Hence: 

< |«i| < hR +1 

. . . 2V3(hR + l) 2 (18) 
3 tii 



6.2.3 Method OVS 

To find the optimal vectors, we define Method OVS. 
Method OVS: 

1. The first step is to search Si the initial set of generator vectors ui and U2- 
Si is the set of vectors having as coordinates the integers (xi,yi) and (^2,2/2) 
verifying System H7I if R > V% and System [TH] if R < V2. 

2. Now, the set S2 should be filtered to keep only reduced and valid vectors. 
Indeed, for each couple of vectors (111,112) in <Si, we should verify: 

2.1. (1*1,142) are reduced, that is they verify System [TO1 

2.2. to check the validity of the coloring, two cases are possible: 
2.2.1 if R > y/2 apply Method VC2. 

2.2.2. otherwise, apply Method VC1. 

3. After the step 2, we obtain the set of valid reduced vectors. Now, the optimal 
vectors are the vectors having the smallest absolute value of their determinant. 

Notice that the search of the optimal vectors can be done by a central unit, that 
distributes the vectors to all nodes. It is also possible that each node in the grid 
computes the two generator vectors. 

We can now evaluate the complexity of VCM that lies in the generator 
vectors computation and in their validity check. 

Theorem 6 Depending on VC method, VCM complexity is in 0(i? 6 ) for Method 
VC1 and 9(i? 4 ) for Method VC2. 

Proof: The vector search phase is in 0(i? 2 ) for each vector. The validity check is 
in 0(i? 2 ) for Method VC1 and 9(1) for Method VC2. ■ 



7 Summary: How to Apply VCM in Practice 

This section summarizes the previous sections. In practice, to apply VCM, we 
start from a set of sensor nodes arranged as a two-dimensional lattice (identi- 
fied by their integer coordinates). In reality, in an actual network, the set of 
neighbors will not be exactly given by the set of nodes within a fixed range R. 
However, notice that a valid /i-hop coloring for a given R, is also a valid h-hop 
coloring for R' < R (although likely non-optimal). Hence, we start by selecting 
the value of R, a radio range such as two nodes at a distance greater than R are 
never neighbors (may be using measurements or neighborhood detection), and 
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a value h. Then, each node proceeds as follows: 

1. Find the optimal valid vectors using the Method OVS. 

2. Each node computes its color by applying either VCM-NCC 1 or VCM- 
NCC 2. 

We can notice that VCM allows each node to know its color in a single round. 

8 Coloring Results with VCM 

Note that further examples of colorings are available externally at |21| . 
8.1 Examples of Vectors 

Table [5] gives for different radio ranges two vectors generating the optimal pe- 
riodic pattern as well as the minimal number of colors obtained by a periodic 
pattern, for both a 2-hop coloring and a 3-hop coloring. The '*' symbol high- 
lights the optimality of the number of colors used. 



Table 2: Vectors generating the optimal number of colors. 



Radio 
range 


2-hop coloring 


3-hop coloring 


vector 1 


vector2 


colors 


vector 1 


vector2 


colors 


1 


(2,1) 


(-1,2) 


5* 


(2,2) 


(-2,2) 


8* 


1.5 


(-3,0) 


(0,3) 


9* 


(4,0) 


(0,4) 


16* 


2 


(3,2) 


(-2,3) 


13* 


(4,3) 


(-3,4) 


25* 


2.5 


(4,3) 


(-1,5) 


23* 


(5,5) 


(-7,2) 


45* 


3 


(5,3) 


(-1,6) 


33* 


(7,5) 


(-8,4) 


68* 


3.5 


(5,4) 


(-6,3) 


39* 


(8,5) 


(-8,5) 


80* 


4 


(7,3) 


(-6,5) 


53* 


(8,8) 


(-11,3) 


112* 


4.5 


(9,2) 


(-6,7) 


75* 


(13,3) 


(-9,10) 


157* 


5 


(9,4) 


(-1,10) 


94* 


(14,4) 


(3,15) 


198* 


5.5 


(9,6) 


(-1,11) 


105* 


(16,0) 


(8,14) 


224* 


6 


(H,4) 


(-9,8) 


124* 


(17,4) 


(-12,13) 


269* 


6.5 


(13,1) 


(-7,11) 


150* 


(-19,0) 


(9,17) 


323* 


7 


(10,9) 


(-4,13) 


166* 


(15,13) 


(-19, 7) 


352* 



We observe that 2-hop coloring of the grid with radio range R = 3 is not equiv- 
alent to 3-hop coloring of a grid with R = 2 in terms of the number f colors (33 
vs. 25). We conclude that the optimal number of colors is not determined only 
by the product h* R, but also the values of h and R separately. 

8.2 Comparison with Other Methods 

Table [3] depicts the simulation results obtained with the VCM for various grids, 
with various radio ranges. The method is compared to a distributed coloring 
algorithm using line/column as priority assignment heuristics. As observed in 
Table [TJ the random priority assignment produces a high number of colors and 
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hence is not presented here. Results are given for 3-hop coloring. For high 
radio range values, the number of nodes should be high enough to allow the 
reproduction of the color pattern. 

Table 3: Number of colors obtained for 3-hop coloring. 



Radio 


Grid size 


Colors 


range 




VCM 


line/ column 


1 


10x10 


8* 


8* 




20x20 


8* 


8* 




30x30 


8* 


8* 


1.5 


10x10 


16* 


16* 




20x20 


16* 


16* 




30x30 


16* 


16* 


2 


10x10 


25* 


30 




20x20 


25* 


33 




30x30 


25* 


33 


3 


20x20 


68* 


80 




30x30 


68* 


83 


3.5 


20x20 


80* 


91 




30x30 


80* 


91 



We observe that VCM provides an optimal 3-hop coloring, for any radio 
range. This is not true for any other priority assignment tested. Moreover, the 
number of colors does not depend on the grid size. 

9 The Vector Method and Real Wireless Net- 
works 

In this paper, VCM has been described for grid topology since this topology is 
used by real applications briefly presented in Section [TJ Notice however, that 
wireless communication may differ from what is expected by the theory that 
often uses simplified models: radio links may be asymmetric, a radio link may 
exist even if the remote node is at a distance higher than the transmission range 
or conversely not exist even if the remote node is in the theoretical radio range. 
That is why, the first step in VCM is to select R such that two nodes that are 
at a distance greater than R are not neighbors. 

Another real aspect in wireless sensor networks is nodes late arrival (because of 
the mobility or in case of late start-up) and nodes disappearance (a node is out 
of battery for instance) . What is the impact of such impairments on VCM? We 
can classify these impairments in two categories: 

1. Radio links disappearance: in this case, VCM always provides a valid col- 
oring. The periodic coloring may still be optimal, as long as the percentage of 
missing radio links is below a given threshold L\. 

2. Radio links appearance: in this case, nodes that should not be neighbors 
(or heard nodes in case of asymmetric links) are. As a consequence, nodes hav- 
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ing the same color may interfere because of these additional radio links. The 
periodic coloring provided by VCM may still be perfectly acceptable by the 
application as long as the percentage of additional radio links is below a given 
threshold L^. 

As a further work, we will evaluate the thresholds L\ and L2 and also study 
how random topologies can be mapped in grid topologies. 

10 Conclusion 

In this research report, we have presented a new method called VCM, the Vector- 
Based Coloring Method, able to provide an optimal periodic h-hop coloring of 
any grid, with h an integer > 1, for any radio range R. This method is easy to 
use: a single round is needed. It suffices to compute the two generator vectors, 
as shown in this paper. Knowing its coordinates within the grid, each node 
deduces its color from a simple computation given in the paper. We have shown 
that this h-hop node coloring is optimal in terms of colors and rounds. We 
determined also an upper and a lower bound for the number of colors needed 
to color an infinite grid. VCM provides the optimal number of colors compared 
to all possible coloring including non periodic ones. Finally, we discussed how 
to apply VCM in real wireless networks. 
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11 Annex 

In this Annex, we group mathematical results and grid properties that are useful 
to study the validity and the optimality of VCM. 

11.1 Gauss Lattice Reduction 

For any pair of vectors u\, ui generating a lattice k{u\,U2), the Gauss lattice 
reduction algorithm provides two reduced vectors v\ , V2 generating exactly the 
same lattice and verifying the System (fTT)|) 



Additional properties are that v\ and v-i are also the two shortest distinct vectors 
of A(tti,M2) and have the same lattice determinant as u\,U2- See for instance 
[T5] for more details. 



11.2 Relation between Number of Hops and Actual Dis- 



Hereafter, we introduce some results related to grid networks. These results can 
be applied to VCM, or any other algorithm. 

Results in this section are inequalities, establishing links between number of 
hops and actual distance. 

Lemma 5 For any point V ofM 2 , there exists a node V of the grid Z 2 such 
that d(V,V) < i/2/2. 

Proof: In the worst case, the node V occupies the center of a grid cell. It is at equal 
distance of two grid nodes that are diagonally opposed. Hence, its distance to 
one of them is equal to V2/2. ■ 

Lemma 6 For any transmission range R > for any grid node U , any node 
V that meets d(U, V) < (R — \/2)h is at most h-hop away from U . 

Proof: We consider the h—1 points of M 2 that allow us to divide the segment [U, V] 
in h equal parts. 

Let Wi be these nodes, with i £ [1, h — 1]. 

For any i £ [l,h — l], let W[ the grid point the closest to Wi. For simplicity 
reason, we denote = U and W' h = V. We have d[U, V) < J2i ~o d(W(, W' i+1 ). 
We have d(W(, W' l+l ) < Wi) + d{W h W l+1 ) + d{W i+1 , W' +1 ). According 

to LemmaEJ we have d(Wi,W{) < V2/2. Hence, we get d{W' i ,Wl +l ) <y/2 + 
d(Wi,W i+ i). By construction, d{Wi,W t+1 ) = d(U,V)/h. 

If d(U,V) < (R- V2)h, then d{W' l ,W' l+1 ) < R. Hence, nodes W t for i e 
{1, 2, . . . , h — 1} constitute a h-hop path between U and V. ■ 

Lemma 7 For any transmission range R, for any two grid nodes U and V , in 
h-hop coloring, if d(U, V) > hR then U and V are at least (h + \)-hop away. 
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Proof: By contradiction assume that, U and V are /i-hop away or less. Let W\ be 
the k — 1 nodes constituting the fc-hop path between U and V, with k < h. Let 
M^i = U, and W), = V. Since nodes Wi are 1-hop neighbors, we have: 

h h 
\UV\ = | J2 WiW i+1 \ < £ d(Wi, Wi+i) < fcfl. 
1=1 1=1 

Hence the contradiction. ■ 

Let U, V be two points of Z 2 and define H(U, V) as the number of hops between 
U and V (it is an integer). For any R > (some inequalities are trivially true 
when R < y/2), the lemma [5] and lemma [7] can be summarized as: 

d(U,V)<(R-V2)h H(U,V)<h (20) 

cf(E7, V) > (/i - l)iZ =^> rl(U,V)>h (21) 

rl(U,V)>h => d(tf,V) > (R-V2)(h- 1) (22) 

H(U,V)<h d{U,V)<hR (23) 

11.3 Bounds on Distance and Number of Hops of Points 
on a Lattice 

Lemma 8 //ui and u 2 are reduced generator vectors of a lattice A(ui, U2), with 
\ux\ < \u2\, then for any vector w such that w = aui + (3u2, and a, /3 € 1 , we 
have \w\ > ja 2 \ui\ 2 , and \w\ > §/3 2 |ui| 2 . 

Proof: Let W the node of coordinates (a,/3). We have: 

\UW\ 2 = a\ Ul \ 2 + /3 2 \u 2 \ 2 + 2aP (ui ■ u 2 ) 

> a \ui\ + ft |t£2 I — 2|a||/3| |ui • «2|- 

Since Mi and W2 are reduced vectors, we can use the property given in the 
System Q32 we get: 

\UW\ 2 > a 2 \ Ul \ 2 +/3 2 \ Ul \ 2 + |a||^||ui| 2 . Hence, 
\UW\ 2 > ((|a|-|/3|) 2 + |a||/3|)| Ul | 2 . 

Notice that this quantity does not change if we change the sign of a or of (5. 
Thus, we assume (a > 0), (0 > 0), and let /(a, 0) = (a - /3) 2 + a/3. 
By a change of variable /? = ^ + A, we get: 

f(a,j3) = |a 2 + A 2 > \a 2 . Similarly, we have f(a,/3) > \fi 2 '. Hence the lemma. 



Lemma 9 Consider any transmission range R > \/2, two reduced generator 
vectors U\ and u-i of the lattice A(u\, U2), and a node W with UW = au\ + j3u2 
for some a and ft in Z. Assume also that the point V\ such that UV\ = U\ is at 
strictly more than h hops from the U . Then: 

if\<A _! M(-R) or \P\ _t M(-R) where /i(i?) = 5^r~7f)j we have: W is strictly more 



than h hops away from U . 
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Proof: Since u\ and U2 are reduced, we can apply Lemma [5] and obtain: 

\UW\ 2 > /W)K| 2 
3 

> — a 2 |ui| 2 , and as a > fJ,(R) 

3 2y/3R . 2| ,2 
" 4 V 3(i?-V2) 



> (_^) 2 | Ul | 2 
i?-V2 



Since the point V\ is strictly more than h-hop away from U , the lemma [5] implies 
by contradiction that |ui| = \UV\\ > (R — y/2)h. It follows that: 

|L%| 2 > i? 2 /i 2 

Applying Lemma we obtain the result. ■ 
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